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ABSTRACT 


A short discussion is given of the reasons for co-rotation of magneto- 
spheres with their planets. It is shown that, if Jupiter's magnetosphere co- 
rotates, the plasma on any magnetic shell must tend to concentrate at those 
parts of the shell which are most distant from the rotational axis, because of 
the large centrifugal force. This would compress the plasma into a discus- 
shaped region about a planetary diameter thick, inclined to the equatorial plane 
by about 7°. The passage of this plasma over the satellite Io offers a new and 
attractive basis on which to explain the influence of Io on the decametric radia- 
tion of Jupiter. It is proposed that the maximum plasma frequency at Io's orbit 
is 40 MHz. Exploration of this possibility leads to a plasma temperature of 
1800°K and a magnetic field of 30 gauss at Jupiter's equatorial surface, with the 
dipole displaced 0.4 planetary radius towards the north rotational pole. Some 
discussion of the difficulties of the theory is given and its future development is 
considered. 
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THE STRUCTURE OF JUPITER'S MAGNETOSPHERE 
AND THE EFFECT OF Io 
ON ITS DECAMETRIC RADIO EMISSION 

1. INTRODUCTION 

The radio emission from the planet Jupiter in the decametric wavelength 
range (< 40 MHz) has been the subject of a great deal of speculation ever since 
its accidental discovery by Burke and Franklin (1955). The demonstration by 
Bigg (1964), that Jupiter's innermost Galilean satellite, Io, has a profound effect 
on the decametric radiation, added a bizarre note to the subject. Several good 
reviews have been published recently, e.g. those by Ellis (1965) and Warwick 
(1967). Most theories depend on the precipitation of electrons from the planet's 
magnetosphere into its ionosphere, but it is hard to see how Io's influence can 
extend over a large enough distance to account for the intensity of the radiation. 

It is observed that the probability of receiving decametric radio noise from 
Jupiter shows two maxima if plotted as a function of the phase angle of Io from 
superior geocentric conjunction, one when the phase is 90° and the other when it 
is 240°. The theory described in this report grew out of an investigation of the 
idea that Io may penetrate the boundary of the Jovian magnetosphere at these 
two positions. 

It is easy to show that this simple theory is untenable. Let us assume that 
the solar wind velocity remains constant as it streams outward to Jupiter's 
orbital distance of 5.2 AU; we shall use 500 km/ sec as a typical value. We as- 
sume also that the particle density falls off in inverse proportion to the square 
of the distance, so that it will be about 1/25 that of the Earth's orbit; selecting 
10 proton-electron pairs per cm 3 for the latter, the density at Jupiter's orbit 

will be cm -3 or approximately 4 x 10 s m“ 3 . The kinetic energy density of 

the solar wind there will then be about 8 x 10“ 11 joule/m 3 . If Io is to penetrate 
the magnetopause at the observed positions of maximum probability of radio re- 
ception, the stagnation point should lie about 5 equatorial planetary radii (5r o ) 
from the center of the planet. Thus, we require the magnetic energy density 
B 2 /2yu. o to be 8 x 10“ 11 joules/m 3 at 5r o . The surface field at the equator of 
the planet is then given by B o = 5 3 B = 125 x 2 n o x 8 x 10“ u = 1.8 x 10 -6 
webers/m 2 or 0.018 Gauss, if we assume a dipole field. The field at 3r o would 
then be 7 x io“ 4 Gauss, which is far too small to contain the high energy elec- 
trons which almost certainly produce the decimetric radiation (Chang and 
Davis, 1962). Estimates based on "reasonable" magnetic field strengths place 
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the magneto spheric stagnation point at 30r o to 70r o (e.g. Field, 1960; Ellis, 1963; 
Carr et al, 1966). 

There is a point which deserves a little attention before we finally reject 
this theory. If the plasma in the magnetosphere co-rotates with the planet, the 
stagnation point in the solar wind is not the point of zero relative velocity. 

Axford (1963) has pointed out that in the Earth's magnetosphere this rotation 
would displace the point of zero relative velocity to the east of the stagnation 
point. In the case of Jupiter, a co-rotating point in the equatorial plane has a 
linear velocity of 500 km/sec at 40r o . Thus there may be no point of zero rela- 
tive velocity at all if some of the earlier estimates of the size of the magneto- 
sphere are correct. This may be important when discussing the shape of Jupiter's 
magneto spheric boundary and circulation inside it. 

2. GRAVITATIONAL AND CENTRIFUGAL FORCES IN 

JUPITER'S MAGNETOSPHERE 

The foregoing discussion leads us to suspect that the rapid rate of rotation 
of Jupiter may be important in determining the distribution of plasma in its 
magnetosphere. This has also been pointed out by Ellis (1965) and Melrose (1967). 
In this section we shall perform some preliminary calculations which confirm 
this suspicion. 

Jupiter rotates extremely fast for so large a body. The period is just less 
than 10 hours, so that the angular velocity is fl = 1.76 x 10~ 4 rad/sec. Since 
the equatorial radius r o is 71,350 km, points on the surface at the equator have 
a linear velocity of 12.5 km/sec due to the rotation. This is almost as large as 
the orbital velocity, 13.1 km/sec. The centripetal acceleration experienced by 
a particle on the surface is thus 2. 2m/ sec 2 . Since the mass of the planet is 
1.90 x 10 27 kg, the gravitational acceleration at the surface on the equator is 
g = 24.9 m/sec 2 , whereas at the pole it is 28.5 m/sec 2 , the polar radius being 
° 14 

66,600 km, or of r o . Thus, even on the surface, at the equator, the centrifu- 
gal force is nearly 1/10 of the gravitational force. This is, of course, the reason 
for the pronounced oblateness of the planet. If the magnetosphere co-rotates, 
the distance r e at which the gravitational and centrifugal forces are equal and 
opposite in the equatorial plane is 2.24 r o . Outside this, the resultant force on 
the plasma, in a coordinate system rotating with the planet, is outward. It is 
obvious that this must have a profound effect on the plasma distribution at, for 
example, Io's orbit, which has a radius of 5.9 r o . We should thus expect Jupi- 
ter's magnetosphere to be different in many respects from that of the Earth, for 
which r is 6.7 Earth radii. 

e 
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3. MAGNETOSPHERIC CO-ROTATION 


It is thus very important to know whether the magnetosphere of the planet 
will co-rotate or not, if we are to discuss intelligently the distribution of parti- 
cles in it. The literature is far from clear on this question and there is still 
some doubt, both from the experimental and theoretical viewpoints, as to whether 
the outer parts of even our own magnetosphere co-rotate or not. No general re- 
view of this subject appears to have been published. Dungey (1958) and Hines 
(1964) have given short accounts of the physical principles involved. 

If we regard the planet as a conducting sphere, of radius r o , rotating with 
angular velocity fi in a magnetic field B , the conduction electrons at the point r 
will experience a Lorentz force -e(fl x V) x B (where -e is the charge on an 
electron). They will redistribute themselves so that the electrostatic field E 
due to the space charge, just cancels this. In the simplest case, where B is 
constant in time at each point, we can then write for the electrostatic potential <t> 
at the point r: 


-E = V4> = (Q x r) x B 


( 1 ) 


The solution of this equation gives, on substituting r = r o , the electrostatic po- 
tential at the surface of the planet. If B varies with time, but in such a way that 
the variation could be attributed to rotation of the source of B with the planet, 
Backus (1956) shows that (1) must be replaced by 

4> = (QxrVA + C (2) 

where A is the (time -varying) vector potential of B and C is an arbitrary 
constant. 

The potential so found at the surface of the planet can be used as the bound- 
ary condition at r = r o for the potential in the magnetosphere, which will be a 
solution of Poisson's equation 


V 2 0) = -— (3) 

^ o 

where p is the charge density and e o the permittivity of free space. By solving 
this equation we could find the electric field E = - V $ at all points in the mag- 
netosphere. The plasma there must then react to this field by undergoing the 
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well-known electric field drift, at a velocity 


v - 


E x B 
B 2 


It is not immediately 


obvious whether this drift velocity will be exactly that required to cause co- 
rotation or not. Nor is it clear how to solve (3), since p depends on the relative 
distribution of positive and negative ions in the magnetosphere, and this in turn 
depends on the field E and so on the rotation. 


Most workers have avoided this difficulty by making assumptions about the 
conductivity of the plasma. The most direct attack on the problem is that of 
Hones and Bergeson (1965). They assume that the plasma is a perfect conductor 
along the lines of force, but has zero conductivity across the magnetic field. 

This implies that there is no component of the electric field along the lines of 
force, i.e. that in the magnetosphere 


E • B = 0 


(4) 


or 



+ 


3A \ 
at ) 


V X A 


0 . 


(5) 


They are able to solve this equation, subject to the boundary condition (2) at the 
surface of the sphere, for the special case in which the magnetic field is that of 
a dipole situated at the center of the sphere, inclined at an arbitrary angle to the 
axis of rotation and rotating with the sphere. The electrostatic potential is 
given by 

$ = ~2 sin 6 [sin 6 cos y - cos 9 sin y cos (<£ - fit)] 


where r, 9 > 4> are the spherical polar coordinates of the point at which <t> is 
evaluated, y is the angle between the rotational and dipole axes, M is the mag- 
netic moment of the dipole and t is the time since the dipole axis was in the 
plane 4> = 0. From (6) Hones and Bergeson are able to show that the average 
drift velocity of a plasma particle is exactly that required for co-rotation, that 
the space charge required to maintain the electric field is negligible compared 
with the particle densities which exist in the Earth's magnetosphere, and that 
the currents due to the particle motions in the field of the rotating inclined 
dipole are not large enough to distort the dipole field appreciably. 
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Presumably, co-rotation begins to fail when the assumptions made in the 
above treatment fail: 

(a) The plasma density may decrease to such a low value, that equation (4) 
no longer holds, so that components of the electric field can exist along 
the magnetic lines of force; 

(b) The plasma may be so dense that it conducts well across the lines of 
force, so that (4) is again invalid and the magnetic field is excluded 
from the plasma; 

(c) The currents due to differential motion of particles of different sign 
may distort the magnetic field appreciably; 

(d) Near the outer boundary, equation (6) will no longer hold, since it will 
have to fit on to the magnetopause and interplanetary potential distribu- 
tion instead of falling steadily to zero at r = oo; 

(e) The plasma density may become too low to maintain the space charge 
required to produce the electric field according to equation (3). 

All these conditions except (b) are likely to be met in the outer parts of 
planetary magnetospheres, but it is very difficult to estimate just where they 
become important. Persson (1966) has suggested that (a) may be true over 
"vast regions” of the magnetosphere of the Earth. It is also well known that 
(c) applies to the Earth's magnetosphere during magnetic storms. (But see 
Hoffman and Bracken, 1967). 

Ferraro has adopted a different approach. In his paper on the rotation of 
the sun (Ferraro, 1937), he showed that, if the plasma does rotate, Maxwell’s 
equations, Ohm's law and the equation of continuity of mass flow combine to 
require that the angular velocity must be constant over the surface traced out 
by the complete revolution of a line of magnetic force about the axis of rotation; 
this is now known as the "law of iso -rotation". It is true, subject to the condi- 
tion that no electric currents flow in the meridianal planes. Recently, Ferraro 
and Bhatia (1967) have extended the theory to planetary magnetospheres, assum- 
ing the conductivity of the plasma to be large but isotropic . They have concluded 
that if there is either no mass flow in the meridianal planes or no azimuthal 
component of the magnetic field, the iso-rotation law holds. They have used the 
mechanical force equation to derive expressions for the pressure (and thus num- 
ber density), but these have only been solved in simple cases. In particular 
there is doubt as to what happens in the region where the centrifugal force ex- 
ceeds that due to gravity. 
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Dungey (1958) has given a useful physical picture of the reason for iso- 
rotation. He points out that, as long as the magnetic field may be considered to 
be "frozen in" to the plasma, any differential rotation of plasma on the same 
line of force must tend to "wind up" the line, so that it might get wound several 
times round the axis. This would be a condition of increasing strain and would 
be resisted by the field. All the material on any line of force must obviously 
rotate with the same angular velocity if the line is to remain in a steady state. 

Presumably, this law breaks down where the plasma is no longer dense 
enough (i.e. a good enough conductor) to hold the magnetic field frozen in. 
Breakdown will also occur if there is appreciable mass flow of plasma along 
the field lines in the meridianal planes, e.g. during the injection of new plasma 
or the ejection of old. Ferraro and Bhatia (1967) themselves suggest that the 
necessary condition for steady rotation is that "the magnetic pressure must ex- 
ceed or be comparable with the centrifugal force.” 

We may note in passing a completely different approach, in which Vlasov 
and Khakimov (1964) have set up and solved the Vlasov equations for the distribu 
tion of electrons and protons in a magnetic dipole field, with rotation of the 
plasma, allowing for deformation of the magnetic field by the currents due to 
particle motions and for the electric field produced by the space charge. Un- 
fortunately, their treatment leads to the result that the ratio of temperature to 
mass of the positive ions must equal that for negative ions. Since this would 
require the ion temperature in the terrestrial magnetosphere to be about 2000 
times the electron temperature, it is clearly unacceptable. Presumably this 
condition results from the absence of a collision term in the Vlasov equation. 

Rotating plasmas have been discussed from the point of view of single par- 
ticle orbit theory by Longmire et al . (1959) and Northrop (1963), among others. 
The discussion by Boyer et al . (1958) of the thermonuclear device "Ixion" con- 
tains several points which are important in understanding magnetospheric co- 
rotation. In "Ixion", a cylinder of plasma is confined in a magnetic field, which 
is uniform except for stronger mirror fields at the ends of the cylinder. A 
radial electric field is applied, which sets the plasma into rapid rotation about 
the axis and hinders escape of particles through the mirror regions by adding 
the centrifugal force to the magnetic forces otherwise opposing escape. The 
reality of the rotation is confirmed by the Doppler shift of spectral lines and by 
the behaviour of the electrical parameters of the device as a whole. The rota- 

E x B 

tional velocity is found to agree with the value calculated for the drift — 

produced by the applied fields. The authors point out, however, that the centri- 
fugal force also produces a drift which is dependent on the mass and charge of 
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the particles, unlike the electric field drift. As a result protons and electrons 
drift with different velocities and there is a net current flowing in the azimuthal 
direction, which causes the plasma as a whole to act as if it were diamagnetic. 
Thus, the magnetic field is reduced in the center of the cylinder; experiments 
whith a probe coil show that it is reduced almost to zero. 

In a planetary magnetosphere we should expect, on the basis of the foregoing, 
that the plasma will rotate primarily because of the drift due to the electric 
field produced by the rotation of the planet in the magnetic field. We should 
expect as a first approximation to find the plasma co-rotating with the planet, 
at least in those regions where its conductivity is high along the lines of mag- 
netic force. Beyond the radius where the gravitational and centrifugal forces 
balance one another, co-rotation, if it takes place at all, will require the mag- 
netic pressure to be at least comparable with the centrifugal pressure. Charge 
and mass -dependent drifts, such as those due to the gravitational, centrifugal, 
pressure and magnetic field gradients, will act to produce a diamagnetic current 
flow in the plasma, which will tend to weaken the magnetic field in the regions 
close to the planet and increase it towards the outside of the region of trapped 
plasma. The net affect will be to distort the lines of force by moving them away 
from the planet, especially in regions where the plasma is most dense (c.f. 
Kendall et al ., (1966); Hoffman and Bracken, (1967)). If plasma accumulates 
to such a density that the centrifugal pressure exceeds the magnetic pressure, 
we may perhaps expect it to break away from the field, possibly in a similar 
manner to that illustrated by Dungey (1958) Figure 5.3. 

4. EQUILIBRIUM OF PLASMA IN A PLANETARY MAGNETOSPHERE 

Let us set up a system of spherical polar coordinates which rotate with the 
planet. A particle of mass m which is at rest in this system at the point r, 

-> r o 

6. -i outside the planet, experiences a gravitational force f = - mg a 

e ° r 2 

and a centrifugal force 

f c - mil 2 r sin 2 6 a r + mQ 2 r sin 6 cos 6 a e 

as shown in Figure 1. These forces may be derived from a potential V, 

f g + f c = - m VV (7) 
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where V = - e — — - Vi Q 2 r 2 sin 2 6 + Constant. (8) 

° r 

Here, a r and a e are unit vectors in the appropriate directions at r, 6 , 4>. 

Positive ions of mass m. 
and charge q i will also ex- 
perience a Lorentz force 
q. (E + v. X B), where E, 

v. and B are all measured 
1 

in the rotating frame . We 
may note that this implies 
that E in the rotating frame 
can have no component per- 
pendicular to the lines of 
magnetic force, for such a 

component would produce a drift of the plasma, with velocity — (Ex B),rela- 

B 2 

tive to this frame, which is contrary to the hypothesis that it is co-rotating with 
the planet. Therefore, in our present frame of reference 

E x B = 0 . (9) 

The ions will move until equilibrium is attained, at which time the net force 
acting is zero. The ion pressure will then be given by the condition 

- Vp. - N. m. VV + N. q. E + N. q. v. x B = 0 (101 

*1 11 1 M 1 M 1 v 

where N i is the number density of ions. Similarly, for electrons 

- Vp - N m VV + NqE + Nqv x B — 0 . (11) 

e e e e M e e 
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We now assume that the plasma consists only of protons and electrons, so 
that q = - q e = e and m g « m. . We also assume that the plasma is electrically 
neutral, so that N. = N g = N. This cannot be strictly true, because the gravita- 
tional and centrifugal fields will tend to separate the protons from the electrons 
and a space charge will be set up to produce an electric field to counteract the 
effect. However, the difference between N g and N. required to maintain this 
field is negligibly small compared with the quantities themselves (see, e.g. 

Alfven and Falthammar, 1963, p. 15) and we may justifiably regard them as 
equal. 

Adding equations (10) and (11) under these conditions, we find 

_ v (Pi + P e ) - Nm. VV + N e (v. - vj x B = 0. (12) 

If we assume a Maxwellian thermal plasma at temperature T, then 

Pi = P e = N k T (13) 

where k is Boltzmann's constant. Noting also that the current density is 

J = Ne(v. - v e ) (14) 

(12) becomes 

- V (2N kT) - Nm ; VV + J x B = 0, 
or if T is constant, dropping the subscript on m. for convenience 

2 k TN V (in N + 


This is merely the equilibrium form of the hydromagnetic equation of motion 
(Alfven and Falthammar, 1963, p. 76), or the momentum equation for a conducting 
fluid (Holt and Haskell, 1965, p. 173) and could have been written down at once 
from either starting point. The derivation given, however, seems more illuminat- 
ing as to the physics of the situation. 
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Strict co-rotation would require v. = v e = 0 in our equations. However, 
as we have already seen, the centrifugal, gravitational and field-gradient drifts 
will be different for protons and electrons, so that there will be a net current. 
This is discussed in Appendix A. 

5. DISTRIBUTION ALONG A LINE OF MAGNETIC FORCE 
If we take the scalar product of equation (15) with B, we get 


V 


|ln N + 


m V \ 
2kT / ' 


B 


0 . 


Thus the gradient of In N + 


m V 
2kT 


is everywhere perpendicular to B, and there- 


fore, along a line of force, 


In N + 


m V 
2kT 


Constant . 


(16) 


Substituting from equation (8) we have 


In N 


mg r 2 

O 

2kT r 


mfi 2 r 2 sin 2 9 


+ 


4kT 


+ Constant . 


(17) 


This is true along a line of magnetic force no matter what the shape of the mag- 
netic field is like. 

6. DISTRIBUTION IN AN ALIGNED DIPOLE FIELD 

We now assume a dipole field, aligned along the rotational axis. Then the 
equation of a line of force is 


r - Lr sin 2 i9 

O 


(18) 
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where L is a constant which gives the distance from the axis, in equatorial 
planetary radii, at which the line of force cuts the equatorial plane. With this 
relation, we can reduce equation (17) to the form 


In N 


mgr 2 
2kT r 


mn 2 r 3 
+ 4kTL r 

O 


+ Constant . 


(19) 


All that remains to be done before we can calculate the plasma number density 
along the line of force is to evaluate the constant. 

An obvious method is to follow the approach of Angerami and Thomas (1965), 
by assuming that the value of N is known at some level in the ionosphere. The 
distances which we are considering are so large that we may consider the iono- 
sphere as lying. at radius r o ; let us assume that the plasma density there is N . 
Substitution in (19) and subtraction then leads at once to the equation 


In N/N 

O 


m io r o 

2kT 



mQ 2 r 2 

O 

4kTL 




( 20 ) 


We assume a temperature of 150°K (Gross and Rasool, 1964) for this preliminary 
calculation. 

It is now a simple matter to substitute the known values and calculate In 
N/N for various values of L. The results forL = 6 are shown in Table 1. 

o 


Table 1 


Plasma distribution along the line of magnetic force L = 6. 


r 

l 

2 

3 

4 

5 

6 

r o 

lo gl0 N ^o 

0 

-144 

i 

-155 

-101 

+17 

+206 
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The conclusion is inescapable, that Jupiter's ionosphere effectively ends long 
before a height of r o above the surface is reached. This is not really surprising, 
for at a temperature of 150°K the scale height of atomic hydrogen is 50 km, so 
r 

that at — = 2 we are 1400 scale heights above the surface. Clearly increasing 

1 o 

the temperature, even by a factor of 100, will not materially affect our conclu- 
sion. We may deduce from Table 1 that the Ange rami -Thomas method is in- 
applicable to Jupiter's magnetosphere, since any plasma there will tend to col- 
lect at the most distant parts of the lines of force and will have no connection 
with the ionospheric plasma. 

This suggests that we should regard the thermal plasma as congregated 
round the furthest point of the line of force from the planet, which is the point 
of minimum potential energy according to equation (8) , being kept there by the 
large centrifugal force and the magnetic field. Let us therefore suppose that 
N = N max at this point, where r = L. We may substitute these values in equation 
(19) and again eliminate the constant by subtraction, to obtain 


In N/N 


max 


m go r o / r o l\ mn 2 r o 2 

2kT \ r h) + 4kTL 



(21) 


It is shown in Appendix B that this corresponds closely to a Gaussian dis- 
tribution of the plasma about the point r = L in the equatorial plane, the width 
of the distribution being proportional to the square root of the temperature T : 


N 


N 

max 


exp 



( 22 ) 


where 


A = 



(23) 
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and s is the distance along the line of force from the equatorial plane. N is 
reduced to 1/e of N v at 

max 


S 



( 24 ) 


Substitution of numerical values, retaining T = 150°K for the time being, gives 
s = 7500 km at L = 6. Thus, the effect of the centrifugal force is to compress 
the plasma into a thin "pancake", only about one tenth of the planetary diameter 
"thick" at L = 6. 

7. RADIAL DISTRIBUTION 


In order to produce a model of Jupiter's magnetosphere on the present 
basis, it is desirable to have some method of estimating N max . Even if we 
could solve the equation of equilibrium perpendicular to the lines of force, 
obtained by taking the vector product of equation (15) with B, we should still 
need a numerical estimate at one point. So far, it has not proved possible to 
solve this equation (cf. Ferraro and Bhatia, 1967). In the absence of a better 
method, we shall take the following simple way out. The magnetic field is the 
agency responsible for holding the plasma in against the centrifugal force, 
through the J x B term. (It should be noted that this implies that the behaviour 
of the plasma with respect to motion across the lines of magnetic force is not 
controlled by collisions, i.e.,the plasma is anisotropic enough to be considered 
"thermal" along the lines of force, but "trapped" against motion across them 
(cf. the "medium density" plasma of Alfven and Falthammar (1963)). We then 
suppose that the field will no longer be able to hold the particles in if the plasma 
density is so great that its rotational kinetic energy per unit volume exceeds the 


mQ 2 


-2 - 


magnetic energy density. The limit thus occurs when l/2 N a 
B 2 

- — , where fj. is the permeability of free space. Introducing the dipole field 
* 


B = B — , we find 
° r 3 


N 


max 



(25) 


This applies in the equatorial plane, where N max occurs and B has its minimum 
value. Equation (25) has been derived by others, e.g. Hines (1964), on a different 
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basis. It has been used to estimate the distance beyond which co-rotation may 
cease, or to a magneto spheric boundary due to rotation, but we shall regard it 
rather differently, as giving the maximum density of plasma which can co- 
rotate at a given radius r. We note that it is an inverse eighth power relation, 
and so will lead to a very rapid decrease of N max with increasing r. 

8. MAGNETOSPHERIC MODEL 1 

In order to calculate the plasma density in our model, we require an esti- 
mate of the magnetic field B o at the planet's equator. We shall arbitrarily 
assume B o = 1 gauss (= 10“ 4 weber/m 2 ), since the effect of changing B o is 
easily visualized by noting from equation (25) that N max is proportional to B o 2 . 
With this figure, and retaining T = 150°K for want of a better estimate, we ar- 
rive at the model shown in Figure 2. This shows a plasma distribution vastly 
different from that surrounding the Earth, and is more reminiscent of Saturn's 
rings than of a magnetosphere. 

9. DISTRIBUTION IN TILTED DIPOLE FIELD 

We now recall from the work of observers of Jupiter's decimetric radiation 
(Morris and Berge, 1962; Roberts and Komesaroff, 1965) that the magnetic 
dipole axis of Jupiter is inclined by an angle of 10° to its rotational axis. On the 
basis of our previous discussion, we would expect the plasma to congregate 
about the points on the lines of force of this tilted dipole at which the minimum 
potential energy occurs. At Io's orbit, the centrifugal term in equation (8) is 
1 10 times the gravitational one; we may therefore get a good idea of the effect 
of the dipole tilt by merely moving the distributions in Figure 2 along the tilted 
lines of force to the points where r sin 6 is at its maximum, i.e. to the points 
furthest from the rotational axis. 

The lines of force of the tilted dipole have the equation 

r = Lr o sin 2 (6 - 10°) 


in the plane in which the dipole axis and the rotational axis lie. (Note that L is 
now the distance in units of r o at which the line cuts the magnetic equatorial 
plane.) Thus, 


r sin 6 - Lr o sin 2 (6 ~ 10°) sin 6 . 
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Differentiation of this with respect to 6 shows that the maximum occurs when 

2 tan 6 - tan (8 ~ 10°) , 

which is quickly solved to give 6 - 96°40'. 

10. MAGNETOSPHERIC MODEL 2 

Thus, the points of maximum plasma 
density will lie approximately in a plane 
tilted 6°40' to the equatorial plane of the 
planet. The resulting model is shown in 
Figure 3, where we have, however, taken 
B to be 30 Gauss in estimating N , 
for a reason which will become evident 
shortly. 

The striking thing about Figure 3 is 
that it shows that, as the plasma disc 
co -rotates with the planet, it will sweep 
across Io twice per relative revolution. 
This offers the first direct clue as to a 
possible explanation of the effect of that 
satellite on the radio emission of the 
planet. 

A prominent feature of the spectro- 
grams of the Io-controlled emission from 
source B is a narrow-band emission 
which rises to a maximum frequency in 
the neighborhood of 40 MHz and then falls 
again (See e.g. Dulk, 1965, Figure 3). 

This offers a striking parallel to the be- 
haviour of the plasma frequency in Io's 
vicinity on our model 2, which would also pass through a maximum as the disc 
swept over Io. This suggests that we should adjust the magnetic field so as to 
make the value of N max at Io's orbital distance appropriate to this frequency. 
Using the well-known relation between electron density and plasma frequency, 
we find N max = 1.24 x 10 4 x (40) 2 = 2.0 x io 7 electrons /cm 3 . Equation (25) 
then requires us to set B o = 30 Gauss, as in Figure 3. 


Figure 2 


JUPITER’S MAGNETOSPHERE - MODEL 2 

DIPOLE TILTED 10° B 0 = 30 GAUSS T = 150° K 



(N in cm 3 ) 


JUPITER'S MAGNETOSPHERE - MODEL 1 

ALIGNED DIPOLE Bq= 1 GAUSS T= 150° K 



(N in cm -3 ) 
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11. PLASMA TEMPERATURE 


With N max adjusted to give the correct maximum frequency, we are in a 
position to estimate the plasma temperature, by matching the time scale of the 
observed frequency change to that expected for the passage of our model distri- 
bution over Io. To do this we require an expression for the distance s of Io 
from the plane of maximum plasma density, i.e. the distance CB in Figure 4. 



Let cp be the angle made by the radius vector to Io with that to the point A at 
which the locus of maximum plasma density crosses its orbit. In the spherical 
triangle ABC, angle, ABC is effectively a right angle, since we want to measure 
s along the line of force; angle BAC is 6° 40'. Then 


whence 


sin 6°40' _ 1 

sin s/r sin cp 


s £ r sin 6°40' sin cp ~ 4.9 x 10 4 sin cp km. (26) 

We may easily show from equation (22) that the plasma frequency varies with s 
as 



40 exp 


( 


- 1.34 x IQ -6 



(27) 
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where the constants are appropriate to Io's orbit, L = 5.9. With (26) this gives: 


f 




40 exp 



sm 


MHz 


(28) 


From the spectra, it takes about 50 minutes for f to change from f max to 0.75 
f max . The plasma rotates with angular velocity 1.76 x IO -4 rad/sec and Io in the 
same direction with angular velocity 0.41 x 10" 4 rad/sec. Thus cp in Figure 4 
increases by 1.34 x 10" 4 rad/sec, or from 0 to 23° in 50 minutes. Equation (28) 

/ sin 2 23° \ 


then gives us 0.75 = exp 


- 3210 


whence we find 


T = 1800° K. 


(29) 


Figure 5 shows the frequency variation computed from equation (28) with T = 
1800°K, compared with that actually observed during four well-defined storms 
from source B (Dulk, 1965). The agreement is good. 

12. MAGNETOSPHERIC MODEL 3. 

Recalculation of the data for Figure 3, using equation (29) and B o = 30 gauss 
leads to the magnetospheric model shown in Figure 6. This has already been 
published in a brief note (Gledhill, 1967). It is worth noting that the temperature 
to which we have been led is comparable with that in the Earth's magnetosphere. 
We should also note that it is unlikely that the plasma density would actually reach 
values of 10 8 or 10 9 cm” 3 , because diffusion across the field and recombination 
would then become important and these have not been considered in our model. 



JUPITER’S MAGNETOSPHERE - MODEL 3 

DIPOLE TILTED 10° B 0 = 30 GAUSS T= 1800° 



N in cm-3 


Figure 5 


Figure 6 
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13. DISPLACED DIPOLE MODEL 


It is obviously a requirement of our model in its present state that the maxi- 
mum frequency of source B should occur when Io passes through the maximum 
density in the plasma plane. Since the magnetic dipole is tilted towards a system 
III longitude of 190° (Roberts and Komesaroff, 1965), this should occur when Io is 
above A m = 100° and 280°. The dependence of the maximum frequency of emis- 
sion on the position of Io in \ m coordinates does not appear to have been stud- 
ied in any detail, though Dulk (1965) has discussed the matter from a somewhat 
different point of view. Examination of a number of published spectrograms of 
the B source in fact shows that the maximum frequency is usually reached when 
Io is within 10° of A m = 230°. Spectrograms of source C also appear to show 
signs of a frequency maximum when Io is near the same value of \ in as for 
source B. Source A spectra are less well-defined in respect of their maximum 
frequency, but a maximum at about A = 150° would be consistent with the 
spectra shown by Dulk (1965, Figure 4). These two positions are equally spaced 
on opposite sides of the plane of the inclined dipole, and are shifted 50° from the 
positions predicted by the model in its present form. 

We may remove this discrepancy, however, if we assume that the dipole is 
not only tilted, but also is displaced from the center of the planet, O, (Figure 7), 
by a distance OD. P and G are the points of intersection of the displaced plasma 
plane with Io's orbit. We require SOP to be 40°, so that, since OP = 5.9 r , 

SO = 4.5 r Q . Thus, since OSD = 6°40', we find OD = SO tan 6°40' = 0.52 r Q . 

A more refined calculation, taking account of the effect of the gravitational field, 
which still acts toward the center of the planet and therefore no longer acts in 
the plane of the disc, and of the changed magnetic field in the disc due to the 
relative tilt, reduces this to 0.4 r Q . This is just at the limit of the displacement 
allowed by the precise observations of Roberts and Ekers (1965). 
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The plasma frequency predicted by this model at Io's orbit is shown in 
Figure 8, together with the frequencies actually observed during seven events 

recorded by the Boulder spectrograph. The 
agreement is reasonably satisfactory for so crude 
a model. If we take into account the observation 
by Roberts and Ekers (1966), that the dipole is 
shifted by 0.07 r o from the planetary axis towards 
= 220°, we should expect a tendency for the 
plasma to concentrate in this direction, and thus 
for the peak at 230° to be higher than that at 150°. 
This is what is observed, as Figure 8 shows. 


Figures 9 and 10 show schematically the re- 
lative configurations of Jupiter, Io, the plasma 
disc and the Earth when the various sources are 
observed to be active. It is seen that in every 
case the radio emission is received at the earth 
when Io is emerging from, or being overtaken by, 
the plasma disc. 

14. COMMENTARY ON THE DISPLACED DIPOLE MODEL 

The main point in favor of the model is that it offers a means by which Io 
can interact directly with the magneto spheric plasma at about the values of 
A- m actually observed. It also predicts the correct type of maximum frequency 
variation, if the radiation is assumed to occur at or near the plasma frequency. 



Figure 8 


CONFIGURATION AT MAXIMUM FREQUENCY 
SOURCE A SOURCE B 


IO'SORBIT 



0 “ 0 “ 
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TYPICAL CONFIGURATIONS 

SOURCE C SOURCE D 



There is no immediate means of explaining the well-known "bursty" nature 
of the radiation observed at decameter wavelengths. The most reasonable basis 
seems to be to assume a "patchy" or turbulent plasma distribution. Alterna- 
tively, one may suppose that Io, by virtue of its gravitational, electric and/or 
magnetic fields, acts by detaching masses from the co-rotating plasma, and that 
these produce radiation as they move away from the planet, through the surround- 
ing plasma. 

The peculiar radiation pattern, by which radio noise is mostly observed 
when Io is around 90° and 240° from superior geocentric conjunction, is perhaps 
explainable on the basis of refraction of the waves by the plasma disc. Some 
very rough ray -tracing shows that there may be a cone above the disc, and one 
below, into which most of the noise is radiated. The cones could reasonably 
make angles of about 6° with the plane of maximum plasma density, more or less 
as required by Gulkis and Carr (1966). Further efforts at ray-tracing do not 
appear to be warranted until the radial distribution of the plasma is put on a 
sounder basis than the present very rough one. 

It is to be expected that the diamagnetic nature of the plasma will result in 
a considerable distortion of the lines of magnetic force inside the disc. One 
would expect them to be bent so as to follow more closely the lines of constant 
plasma density in Figure 8. Some tendency to bend backwards, due to the ten- 
dency of outward moving plasma to conserve its angular momentum, is also to 
be anticipated. As a result, it appears to be impossible at present to make any 
valid predictions about the polarization and Faraday rotation produced by the 
model. 
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The weakest features of the model are undoubtedly the radial distribution 
and the relatively high plasma density required. We shall examine each of these 
a little more closely. 

(a) The radial distribution 

If the plasma were assumed to be perfectly conducting isotropic thermal 
plasma, it would exclude the magnetic field from its interior and there would 
be no force to hold it in against the centrifugal one. The model depends on 
the assumption that motion along the lines of force is sufficiently easy for 
thermal equilibrium to be maintained by collisions and the electric field to 
be removed by the high conductivity, while motion across the lines of force 
is sufficiently difficult to allow the plasma to be regarded as a non-conductor 
in that direction. Calculation on the basis of the model gives the electron- 
ion collision frequency at Io's orbit as 6.5 kc/sec, while the electron gyro- 
frequency is 65 kc/sec; an electron completes, on the average, 10 orbits 
round a line of force before suffering a collision. The plasma is thus of 
"medium density," in the sense of Alfv£n and Falthammar (1963). This is 
marginal for stability and would lead to fairly rapid diffusion across the 
field lines. This in turn would require a suitably strong source of replace- 
ment plasma to retain the high density required. It may also be noted that 
the electron gyrofrequency, o> e , being proportional to B, decreases out- 
wards as r -3 , while the electron-ion collision frequency, v ei , is roughly 
proportional to N and so varies as r -8 . Thus the ratio v ei /w e varies as 
r“ 5 and from this point of view the model distribution is reasonably stable 
outside Io's orbit, but would certainly not be so inside 4r o , where v e . > 
a) e /2rr , i.e. the plasma would take on "high density" characteristics. 

The force holding the plasma is in fact the J x b force in equation (15) 
and the current density J is due to mass- and charge -dependent drifts as 
discussed in Appendix A. This azimuthal drift current must produce a mag- 
netic field which opposes the applied dipole field inside the current ring and 
this effect will be most pronounced in the denser plasma. This is the field 
which produces the distortion of the lines of force already mentioned in 
Section 14. Rough calculations show that it would reduce the applied field 
very seriously inside 5r o , so that the model cannot be even approximately 
correct inside this distance. 

Finally we may note that our assumption of a constant temperature is 
not necessary in view of the assumed difficulty of radial motion of electrons 
across the lines of force. 
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Quite obviously the problem of the radial distribution is a complex one, 
and demands serious study before the model can be regarded as even ap- 
proximating to reality. 

(b) The high plasma density 

The density required at Io's orbit by the model is 2 x 10 7 electrons and 
the same number of protons per cm 3 . 

If we believe the model as far in as 5r o this rises to 8 x 10 7 . Such 
high densities are not encountered in the Earth's magnetosphere and one 
instinctively shrinks from accepting them for Jupiter's. On the other hand 
we have no first-hand knowledge of the Jovian magnetosphere and, given 
the larger magnetic field and an adequate source, say from the solar wind, 
there seems to be no fundamental reason why high densities should not 
exist there. 

Calculation suggests that a density of the order of 2 x io 7 cm -3 over 
the , volume involved should give out a detectable amount of radiation in the 
H a line, though the detection may be difficult. It also appears that the radia- 
tion produced by free -free transitions may be detectable in the 6-20 cm 
wavelength range, though here again an examination of published papers 
suggests that it would be difficult to distinguish such radiation from the 
thermal component coming from the planet itself. 

(c) High energy plasma 

It has been assumed in the theory that only thermal energies are in- 
volved. The microwave observations show that there are electrons with 
energies of several MeV in the region L = 1.5 - 3 and it would be surpris- 
ing if there were no particles further out with energies in excess of the 
thermal ones. Such particles, mirroring in the normal way but with their 
orbits modified by the centrifugal field, might offer a source for the radia- 
tion observed to come from the region of L = 6.5 by McAdam (1966). They 
would also serve to maintain conduction along the lines of force above and 
below the plasma disc, which may be necessary for its stability and co- 
rotation. 

(d) Radiation mechanisms 

No serious attention has yet been given to this aspect of the problem. 

It is well known that a plasma does not, in general, radiate electromagnetic 
waves at the plasma frequency. With respect to the Io-controlled emission, 
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however, it may be significant that the streaming velocity of the co-rotating 
plasma past Io, 55 km/ sec, is 11 times the velocity of sound, 5 km/ sec; it 
is lower than the Alfv<5n velocity, which is 72 km/sec at the maximum of 
the plasma density at Io's orbit and is greater than this elsewhere at L = 5.9. 
Under these circumstances a modified sonic shock wave would be expected 
to precede Io in its orbit (since the plasma rotates faster than Io does). A 
mechanism such as that proposed by Tidman (1965) could then be invoked to 
account for the electromagnetic radiation. It remains to be seen whether it 
would be capable of explaining the observed intensity of the Jovian radiation. 

That part of the radiation which is not controlled by Io may similarly be 
attributed to turbulence at, or the breaking away of clouds of plasma from, 
the boundary region between the co-rotating and the non-co-rotating plasma. 
This must, on the model, occur outside Io's orbit, and therefore, at a lower 
plasma frequency than that in Io's vicinity. Gruber (1967) has recently shown 
that, whereas Io exercises almost complete control at frequencies above 30 
MHz, there is little evidence for such control below 25 MHz, except for a 
short range at 18-19 MHz. It may be significant that the higher range of 
Io-control centers on 36 MHz, which is the second harmonic of the 18 MHz 
band, as Gruber points out. This would seem to be consistent with Tidman' s 
mechanism, which produces radiation at both the plasma frequency and its 
second harmonic. Much more work is necessary to examine this possibility. 
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Appendix A. 

Particle drifts in Jupiter's magnetosphere 


The drift velocity (normal to B) of the guiding center, produced by a force 
F, is (Alfven and Falthammar, 1963, Ch. 2) 


1 _ 

v - F x B 

qB 2 


(Al) 


where q is the charge on the particle. There will thus be the following drifts 
imposed on the ions in the Jovian magnetosphere: 

(a) The pressure-gradient drift, 


pi 


1 1 _ 

;B 2lT 


VP. x B 


kT 

eB 2 


V (In N) x B. 


(A2) 


(b) The gravitational drift 


m. 


v 


eB 2 ° r 2 r 


a x B . 


(A3) 


(c) The centrifugal drift 


1 1 


Vci r>2 2 m 


eB 2 


. fi. 2 V ( r 2 sin 2 6 ) x B. 


(A4) 


(d) The magnetic field drift, produced by the gradient and curvature of 
the field lines: 


Bi 


1 m i 

eB 2 B 




VB x B 


(A 5) 
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where v t and v |( are the components of v. normal and parallel to B 
respectively. 

(c) The electric field drift 


1 -. 

— E x B 
B 


(A 6) 


where the subscript " i" is unnecessary since the velocity is independent 
of the charge and mass. There are corresponding expressions for the 
electrons. 

It is instructive to evaluate these drift velocities for our model, at 
Io's orbit. The results are given in Table Al. All drift velocities are 
in the azimuthal direction. A positive sign indicates the direction of 
co-rotation and B has been taken as directed southwards at Io's orbit. 

In calculating the magnetic field drift, the second term, due to curva- 
ture of the lines of force, has been neglected. This is acceptable for 
the undistorted dipole field but may be unrealistic if the distortion is 
taken into account. The electric field used in computing the drift (e) 
was taken as that found by Hones and Bergeson (1965). 


Table Al 


Drift 

Proportional 

to 

Metres/sec 

Ions 

Electrons ’ 

(a) Pressure gradient 

r 2 

+ 2.2 x 10" 4 

- 2.2 x 10"4 

(b) Gravitational 

r 

- 5.1 x io~ 4 

+ 2.8 x 10" 7 

(c) Centrifugal 

r 4 

+ 9.7 x l(T 3 

- 5.0 x 10 -6 

(d) Magnetic 

r 4 

+ 1.5 x io~ 3 

- 8.0 x 10“ 7 

(e) Electric 

r 

+ 7.3 x io 4 

+ 7.3 x 10 4 


This illustrates clearly the overwhelming importance of the electric field 
in producing co-rotation and shows how closely one may approximate by regarding 
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the electrons and ions as both co-rotating with the same velocity. The net drift 
velocity of the ions in the co-rotating frame is about 10 _2 m/sec and that of the 
electrons is 2 x 10” 4 m/sec. Thus, most of the current is due to the ion drift. 


Appendix B. 

Gaussian distribution of plasma 


The diagram shows the line 
of force L, and the radius vector 
r to the point P. Provided that 
the latitude k is small, we may 
write for the distance s along 
the line of force from Q to p. 



s £ Lr o k . 


(Bl) 


The equation of the line of force is 


r - Lr sin 2 6 - Lr cos 2 k 

O O 


% Lr o (1 - \ 2 ) 


so that 



'V/ 


L (1 - ^ 2 ) 




L 



(B2) ' 


using (Bl). 
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Thus, in equation (21), we may write, since 




(B3) 


and 



(B4) 


(21) then becomes 


In N/N 


max 


a; 

Ay 


mg oS 2 

2kTL 3 r 

O 


mQ 2 3 s 2 
4kTL 3 


Ay 

Ay 


m 

2kT 


gp 

r o L3 



(B5) 


Then we can write 



(B6) 


where 



(B7) 


which are (22) and (23). 
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